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3 359,409 

CORRELATION FUNCTION APPARATUS 
Hugh L. Dryden, Deputy Administrator of the National 
Aeronautics and Space Administration, with respect to 
an invention of Richard M. Goldstein, Los Angeles, 
Calif. 

Filed Apr. 9, 1963, Ser. No. 271,824 
7 Claims. (CL 235—181) 

This invention relates to an improved method and ap- 
paratus for developing correlation functions of signals. 

The invention described herein was made in the per- 
formance of work under a NASA contract and is subject 
to the provisions of Section 305 of the National Aero- 
nautics and Space Act of 1958, Public Law 85-568 (72 
Stat. 435; 42 USC 2457). 

The development of correlation functions is of signifi- 
cant interest in the field of communications where, for 
example, it is desired to identify a signal contained in a 
very noisy background, as by obtaining the power spec- 
trum or signal power density versus frequency, of the sig- 
nal. A well known technique for determining the power 
spectrum of such a signal is to subtract the power spec- 
trum of the noise plus the signal from the power spectrum 
of the noise taken alone. The difference in the power 
spectra constitutes the power spectrum of the desired 
signal. 

A well known technique for developing the power 
spectrum of a signal consists of first determining the auto- 
correlation function of the signal and then finding the 
Fourier transform of the autocorrelation function. This 
technique is discussed in some detail on page 7 of the book 
entitled “The Measurement of Power Spectra/' Blackman 
and Tukey, Dover Publications, 1958. 

The correlation function with which this invention is 
concerned is identical to that discussed in U.S. Patent 
2.854,191 and is defined in general by the following equa- 
tion: 
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In the foregoing equation the terms (MO and f 2 (t ) may 
be different function^ of time, in which case the equation 
defines the cross correlation between the functions. On the 
other hand, the terms may refer to the same function, in 
which case the foregoing equation reduces to 
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which defines the so-called autocorrelation of the func- 
tion / 2 (7). In either case, the left-hand side of the equa- 
tion is a function of r only since / is eliminated by the 
process of integration. 

The foregoing equations define the correlation only for 
the limiting in which T, the period over which the inte- 
gration is taken, approaches infinity. In practice, how- 
ever. interesting and valuable results are obtained in 
which T, while large compared with r, is finite. That is to 
say, it is of advantage to define a “short term” cross cor- 
relation R xl2 (Lr) and a “short term” autocorrelation 
R x11 U,t) respectively as 
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Evidently the quantity R x i 2 (M) is a function of the 
10 time t, as well as the lag r, and so varies as time pro- 
gresses. However, at any instant t, it is a good approxi- 
mation to Q x12 (t), provided only that T»r. To the same 
extent, and within the same restriction, R x n(Lr) is a good 
15 approximation to Q xll (r). 

In the Equations 1 and 2, the period of integration 
commences in the past and extends into the future. It is of 
course impossible to instrument this entire integration 
process with physical means. Therefore, in the Equations 
20 3 and 4, the integration is stated to commence in the past 
and to extend to the present. The integrals as thus changed 
are now physically realizable, though the change in form 
involves no change in their values. 

The computation of correlation is of chief interest in 
25 the field of communications when the correlation (auto- 
correlation or cross correlation as the case may be) for 
one value of the lag r is compared with that for other 
values of the lag r. It is desirable to effect such a com- 
3Q parison continuously and on a running basis as the sig- 
nals whose correlation is being determined progress and 
it therefore becomes necessary to compute the correlation 
for each of a large number of values of the lag r with 
great rapidity. The result of the computation may then be 
35 plotted as a running graphical function of time which 
serves to reveal features of the signals MO and / 2 (r) 
which are otherwise obscured. 

As a practical matter, according to conventional tech- 
niques, the correlation function is developed by measur- 
40 ing signal amplitudes at small finite intervals (r), e.g. 
15 microseconds, over a relatively long time period (T), 
e.g. 1 hour. Then, a first plurality of products is devel- 
oped, each of whose factors are signal amplitudes spaced 
45 in time by the interval r and the average of the:e prod- 
° ucts is determined. Then a second plurality of products 
is developed each of whose factors are signal amplitudes 
spaced in time by the interval 2 r and the average of these 
products is determined. Similarly, the averages of prod- 
50 ucts whose factors are signal amplitudes spaced in time 
by the intervals 3r, 4r, . . . tit, where n may, for ex- 
ample, equal 44 are developed. The relationship between 
these developed averages products and the correspond- 
ing intervals or time delays is of course the correlation 
55 function. Where the factors or signal amplitudes all re- 
late to the same time varying signal, the function is called 
the autocorrelation function, and is principally to this 
function that the discussion herein will be directed, al- 
though it is to be understood that the invention is equally 
applicable to the development of a cross correlation func- 
tion. 

In order to insure that the developed autocorrelation 
function depends only on the recurring characteristics of 
65 the signal and not on its random variations, the exemplary 
quantities of r~ 15 microseconds and T=1 hour rep- 
resent practical choices. From the nature of the calcula- 
tions required to develop the autocorrelation function, 
it should be apparent that the use of a high speed digital 
70 computer is called for. However, because of the exceed- 
ingly great number of calculations which must be per- 
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formed to develop the function for a signal period as long 
as one hour, a very great amount of computer time, pos- 
sibly on the order of 1,000 hours is required. 

In view of the amount of computer time required to 
develop a correlation function by conventional tech- 
niques, it is an object of the present invention to provide 
an improved method and apparatus for continuously 
developing a correlation function as the signal or signals 
are received so as to completely develop the function 
during this receiving period. 

It is a further object of the present invention to pro- 
vide an improved method and apparatus for developing 
a correlation function which requires considerably less 
complex circuitry than is required by apparatus operat- 
ing in accordance with prior art techniques. 

Briefly, the invention herein is based on the recogni- 
tion that a correlation function can be continuously 
developed by considering, at an appropriate sampling 
rate, only the algebraic sign of the signal or signals, 
rather than the amplitude thereof and that by sequential- 
ly applying binary signals representative of the sign to 
delay means consisting of many multiples of basic de- 
lays, many such binary signals, each with a different de- 
lay are simultaneously made available to thereby permit 
products of signals separated by a certain time delay to 
be continuously taken. 

In a preferred embodiment of the invention, a signal 
whose autocorrelation function is desired, is applied to 
a limiter which operates to develop binary signals (rep- 
resentative of + 1 and — 1 ) in accordance with the alge- 
braic sign of the signal. The binary signals are sampled 
at an appropriate rate (Nyquist rate) and are successively 
applied to a delay line consisting of a plurality of identi- 
cal delay stages. The delay line is provided with a plu- 
rality of output taps each of which is connected to the 
line between a pair of adjacent delay stages. A different 
logical comparer is connected to each of the taps and to 
the input terminal of the delay line and functions to find 
the product of the applied binary signals, which it should 
be recognized is equivalent to comparing the applied 
undelayed binary signal with the applied binary signal de- 
layed by a specific amount of time for identity. For every 
positive comparison, i.e. where the binary signals applied 
to the comparer are identical, an accumulator connected 
to the particular comparer is incremented by one. By 
counting the total number of signals applied to the delay 
line and subtracting this count from twice the count in 
the accumulator, the numerical difference between the 
number of positive and negative comparisons and con- 
sequently the average of the products of the signs of the 
signal for the specific delay time can be determined. The 
plurality of averages provided by the comparers and ac- 
cumulators associated with the respective taps constitute 
a distorted autocorrelation function due to the distortion 
introduced by the limiter. The true autocorrelation func- 
tion can be represented by R xU which is equal to 

sin B» u ) 

where R yll represents the distorted autocorrelation func- 
tion. The relationship between the true and distorted 
autocorrelation functions is discussed in the text entitled, 
“Random Signals and Noise/* by Davenport and Root, 
McGraw-Hill, 1958, on page 107. 

The novel features that are considered characteristic of 
this invention are set forth with particularity in the ap- 
pended claims. The invention itself both as to its organ- 
ization and method of operation, as well as additional 
objects and advantages thereof, will best be understood 
from the following description when read in connection 
with the accompanying drawings, in which: 

FIGURE 1 is a block diagram of a typical prior art 
autocorrelation system showing the waveforms and func- 
tions obtained at each step in the development of the auto- 
correlation function; 
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FIGURE 2 is a block diagram of the present invention 
showing the waveforms and functions obtained at each 
step in the development of the autocorrelation function; 
and 

5 FIGURE 3 is a schematic diagram illustrating a pre- 
ferred embodiment of the present invention. 

Attention is initially called to FIGURE 1 which repre- 
sents a typical prior art system for compiling the power 
spectrum of a signal. The incoming signal 10, which can 
comprise a noisy signal with or without the signal whose 
power spectrum is ultimately desired, is applied to an 
autocorrelator 12 to obtain the autocorrelation function 
14 of the signal. The autocorrelation function is obtained 
by averaging the products of signal amplitudes spaced in 
15 time by an integral number of finite time delays (r). 
That is, assuming the quantities P 0 , Pi, P2, * • • Pn rep- 
resent the amplitude of the signal 10 at times 0, r, 2 r 
... nr, the average of the products P 0 Pi, PiP 2 , P 2 Ps 
. . . Pn-iPn represents the initial point on the autocor- 
20 relation function 14. Similarly the average of the products 
P 0 P 2 , F1P3, P2P4 • • • Pn— aPn represents a second point, 
etc. The autocorrelator 12 performs the necessary calcu- 
lations to obtain all points of the autocorrelation function 
14. In order to develop the power spectrum of the signal 
25 from the autocorrelation function 14, the autocorrelation 
function is applied to a calculator 16 which in turn yields 
the Fourier transform 18 of the auto-correlation function 
14. As is well known, the Fourier transform of the auto- 
correlation function constitutes the power spectrum, i.e. 
30 power density versus frequency, of the signal 10. 

Attention is now called to FIGURE 2 which illus- 
strates, in block diagram form, the technique for develop- 
ing the autocorrelation function and ultimately the power 
spectrum in accordance with the invention. Essentially, 
35 the technique according to the invention differs from the 
prior art technique in that instead of finding the products 
of the amplitude of the signal at points spaced in time, the 
products of the algebraic signs of the signal are utilized. 
More particularly, again assume it is desired to develop 
40 the power spectrum of an arbitrary signal 10. The signal 
10 is applied to sampler means 20 which in turn func- 
tions only to provide the algebraic signs of the signal 10 
at an appropriate rate. More particularly, note the output 
signal 22 derived from sampler means 20. Gutput signal 
22 consists of a series of equal amplitude pulses which are 
0 either positive or negative ( -f-1 or —1) in accordance 
with the algebraic sign of the signal 10. Output signal 22 
is in turn applied to a delay circuit 24 which is provided 
with a plurality of output lines 26 each representative of 
50 a different integral number of basic delay periods. That is, 
if the signal 10 was being sampled at a 15 microsecond 
(r) rate, no time delay would be associated with a first 
output line of the means 24, a time delay equal to r would 
be associated with a second output line, a time delay equal 
55 to 2 r would be associated with a third output line, etc. 
The signal patterns respectively appearing on several of 
the delay circuit output lines 26 are generally identified 
by the designating numeral 28. The output lines 26 of 
the delay circuit 24 are in turn respectively connected to 
60 a plurality of comparer circuits 30 each of which func- 
tions to obtain the product of the signal appearing on the 
first output line and the signal appearing on the unique 
output line to which it is connected every 15 microsec- 
onds. It should be apparent that the product will be ei- 
65 ther +1 or — I. The comparer circuits 30 are in turn con- 
nected to an averager 32 which functions to continually 
maintain an average of all of the products respectively 
associated with each of the various delay times. More 
7 particularly, the averager 32 yields the average of all 
products having factors spaced in time by 15 microsec- 
onds, the average of all products having factors spaced in 
time by 30 microseconds, etc. These averages are con- 
tinuously taken so that after a period of, for example one 
75 hour, ail points of a distorted autocorrelation function 
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Ryu, represented by the curve 34, are known. The dis- 
tortion is of course due to the fact that only the algebraic 
signs of the signal 10 were considered rather than the 
exact amplitudes of the signal 10. However, it has been 
found that the distorted signal R yl i can be modified to 
arrive at the true autocorrelation function R^i by finding 



This calculation can be performed by calculator 36 to 
thereby provide the true autocorrelation function Rjai, 
represented by the curve 14 at the output thereof. The 
true autocorrelation R 2 n can then be applied to a calcu- 
lator 16 to develop the Fourier transform IS according 
to conventional techniques. Calculator 36 can, for ex- 
ample, comprise a special purpose circuit or any conven- 
tional general purpose digital computer. 

From the foregoing, it should be apparent that the 
technique described by FIGURE 2 differs substantially 
from the prior art technique described by FIGURE 1 
in that the autocorrelation function is compiled concur- 
rent with the reception of the signal 10 so that at the 
end of the sampling period of the signal IS, the com- 
pleted autocorrelation function is available. The con- 
tinual real time compilation of the autocorrelation func- 
tion avoids the necessity of devoting a substantial amount 
of computer time to its compilation. The reason, of 
course, which permits the continual compilation of the 
autocorrelation function, as shown in FIGURE 2, is be- 
cause only the algebraic signs of the signal 10 are con- 
sidered rather than the exact amplitudes of these signals 
as shown in FIGURE 1 and in addition because the delay 
circuit 24 simultaneously provides signals having many 
multiples of a basic delay. In addition to eliminating the 
requirement of utilizing a substantial amount of com- 
puter time, it should be apparent that considerably less 
complex circuitry is required to perform the operations 
necessitated in FIGURE 2 than is required to perform 
those essential to FIGURE 1. 

Attention is now called to FIGURE 3 which illustrates 
in schematic form a preferred embodiment of the inven- 
tion illustrated in block form in FIGURE 2. More par- 
ticularly, the sampler means 20 is provided with a limiter 
circuit including a resistor 4© connected to ground 
through a pair of oppositely poled paralleled diodes 42 
and 44. The signal 10 whose autocorrelation function is 
desired is applied to input terminal 46 connected to the 
resistor 40. The input terminal 48 of a flip-flop 50 is 
connected to the junction 41 defined between the resistor 
40 and diodes 42 and 44. The flip-flop 50 is provided with 
true and false output terminals which are respectively 
connected to the inputs of And gates 52 and 54. The sec- 
ond input to each of And gates 52 and 54 is connected to 
the output of a conventional clock pulse generator S6 
which functions to provide clock pulse signals at a 
specified rate (Nyquist rate, herein assumed to be every 
15 microseconds). 

The outputs of the And gates 52 and 54 are respec- 
tively connected to delay lines 58 and 60. Each of the 
delay lines 58 and 60 includes a chain of delay elements, 
each element capable of delaying a signal applied to its 
respective delay line by an identical finite period equal to 
r (15 microseconds). Taps are connected along the length 
of each delay line to permit ^application of the signals 
applied to the delay line to appropriate ones of the com- 
parer circuits 30. The nomenclature utilized to identify 
the various taps is as follows: the taps connected to the 
delay line 58, which in turn is connected to the true out- 
put terminal of the flip-flop 50, will be identified by the 
upper case letter “A” together with a subscript indica- 
tive of the number of delay stages between the identified 
tap and the input to the delay line. Consequently, the 
initial tap connected to delay line 58 will be identified 
as Aq meaning of course that an undelayed signal can be 
sample at this tap. Similarly, the next tap connected to 
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delay line 58 will be identified as Ai meaning that the 
signal sampled at this tap will be delayed by one finite 
period equal to r. The next tap will be identified as A s 
meaning that a signal sampled at this tap will be delayed 
g two finite periods equal to 2r. The taps connected to de- 
lay line 6 © will be similarly identified except however 
inasmuch as this delay line is connected to the false out- 
put terminal of flip-flop 50, the symbol A will be utilized 
in lieu of the symbol A. 

IQ In order to compile the average of the products for 
each specific time delay, a plurality of comparer circuits 
30 is provided. The individual comparer circuits 30 will 
be identified by the numeral 30 together with a subscript 
0 , !, 2 , . . . n to indicate the integral number of finite 
15 periods constituting the time delay between the binary 
signals applied to it. That is, comparer circuit 30 3 for 
example is dedicated to comparing or finding the product 
of all signals separated by an interval equal to three finite 
periods (3r). Each of the comparer circuits 30 includes 
20 a pair of And gates 62 and 64 whose outputs are re- 
spectively connected to the inputs of an Or gate 66 . The 
output of clock pulse generator 56 is connected to the 
input of each of the And gates 62 and 64 of each of the 
comparer circuits 30. Similarly, the delay line taps A<> 
25 and Ao are respectively connected to the inputs of the And 
gates 62 and 64 of all of the comparer circuits 30. Addi- 
tionally, the delay line taps A x and A x are respectively 
connected to the And gates 62 and 64 of comparer circuit 
30i, the taps A 2 and A 2 to the inputs of And gates 62 
and 64 respectively of the comparer circuit 3© 2 , etc. 

The outputs of the Or gates 66 of each of the comparer 
circuits 30 is respectively connected to an accumulator 
68 forming part of the averager 32. The accumulators 
68 perform a counting function. That is, each of the ac- 
35 cumulators 68 functions to count the number of out- 
put signals generated by the Or gate 66 connected to the 
input thereof. The accumulator 68 0 is connected to the 
comparer 3 © 0 and its output is connected to the input of a 
plurality of arithmetic circuits 20 , also forming part of 
4U averager 32, each of which, is in addition connected to 
the output of a different unique one of the accumulators 
68 . The output of each of the arithmetic circuits 70 rep- 
resents a point on the distorted autocorrelation function 
curve 34 previously discussed in connection with FIG- 
45 URE 2. By applying the output from each of the arith- 
metic circuits 70 to previously mentioned calculators 36 
to perform the required correction, the true autocorrela- 
tion function R kU can be provided. The output of the 
correction circuits 36 is in turn applied to calculating 
means 16 which functions to find the Fourier transform 
of the autocorrelation function R xll which, as was previ- 
ously noted, results in the development of the power 
spectrum of signal li. 

In the operation of the embodiment of FIGURE 3, the 
55 signal 10 is applied to the input terminal 46 of the 
sampler means 20. So long as the signal 10 is positive, a 
current will exist in resistor 40 and diode 42 and accord- 
ingly establish a positive voltage at the junction 41 equal 
to the drop across diode 42. On the other hand, if the 
60 signal 10 is negative, a negative voltage of equal mag- 
nitude will be established at junction 41 by virtue of the 
drop across diode 44. By connecting junction 41 to the 
input terminal 48 of flip-flop 50, the flip-flop is caused 
to assume one state, for example the true state, in re- 
85 sponse to a positive voltage level at junction 41 and its 
other or false state in response to a negative voltage level 
at junction 41. When the flip-flop 58 is in its true state, 
an output signal will be provided by And gate 32 in re- 
sponse to the application of a pulse thereto by clock pulse 
generator 56. On the other hand, when flip-flop 50 is in 
its false state, the And gate 54 will provide an output 
signal in response to the application of a pulse thereto 
by the clock pulse generator 56. The pulses so generated 
75 by the And gates 52 and 54 will respectively travel down 
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the delay lines SB and 6®, being delayed by a finite 
amount by each delay element. 

One of the And gates 62 or 64 of the comparer circuit 
30 0 will, of course, provide an output pulse each time a 
clock pulse is generated inasmuch as a pulse must concur- 
rently appear at one of the taps A 0 or A 0 . Consequently, 
Or gate 66 will provide an output signal which will in- 
crement the accumulator 68 0 by one. In other words, the 
accumulator 68 0 will maintain a running count of the 
total number of pulses applied to the delay lines 58 and 
60 . On the other hand, comparer circuit 3®! will provide 
an output signal to accumulator 68i only when the product 
of two signals spaced in time by a time delay equal to r 
is one. More particularly, by connecting taps A 0 and Ai 
to the inputs of And gate 62 of comparer circuit 30x and 
tap A 0 and A x to the inputs of And gate 64 of comparer 
circuit 3®i, an output signal will be provided by Or gate 
66 of the comparer circuit only when the product of the 
signals spaced in time by r is equal to -f 1. The ac- 
cumulator 68! therefore keeps a running count of the 
number of identical comparisons between signals spaced 
in time by a time delay equal to r. By knowing the total 
number of signals applied to the delay lines 58 and 60 
and by knowing the count in each of the accumulators 
68 1 ~68 n the number of products equal to +1 and —1 
and consequently the. average of the products of signals 
spaced in time by each integral number of finite periods r 
can be determined. The arithmetic functions necessary to 
arrive at these averages are performed in the arithmetic 
circuits 70. As a consequence, the outputs of the various 
arithmetic circuits 70 represent the average product for 
the particular time delay with which the arithmetic cir- 
cuit is associated. As pointed out, by operating upon these 
average products in calculator circuits 34 and by then 
finding the Fourier transform of the resulting autocor- 
relation function, the desired power spectrum of the 
signal 10 is developed. 

From the foregoing, it should be apparent that a method 
and an apparatus for more easily compiling the correla- 
tion function, be it the autocorrelation or cross correla- 
tion function of a signal (or signals) has been disclosed. 
More particularly, the specification has herein taught a 
method and apparatus for continuously compiling a cor- 
relation function so that the correlation function is im- 
mediately available at the termination of the sampling 
period over which the compilation of the correlation 
function is desired. The continuous compilation of a cor- 
relation function concurrent with the sampling of the 
signal for which the function is desired of course re- 
presents a considerable improvement over the previously 
mentioned prior art techniques which necessitate con- 
siderable computing time to arrive at the same results 
arrived at by utilization of the invention. 

I claim: 

1. Apparatus for generating the autocorrelation func- 
tion of a specified signal comprising; first means for suc- 
cessively developing binary signals, at a specified sampling 
rate, in accordance with the algebraic sign of said specified 
signal; a plurality of second means for respectively com- 
paring binary signals representative of the sign of said 
specified signal spaced in time by predetermined respec- 
tive amounts for determining a distorted autocorrelation 
function of said specified signal; and means for correct- 
ing said distorted autocorrelation function. 

2. Apparatus for generating the autocorrelation func- 
tion of a specified signal comprising first means for 
developing binary signals in accordance with the algebraic 
sign of said specified signal; a plurality of delay means 
each capable of delaying a signal applied thereto by a 
different specified amount; sampling means for applying 
said developed binary signals to all of said delay means 
at a specified rate; a plurality of comparing means each 
of which is connected to a different unique one of said 
delay means for continuously comparing for identity the 
delayed binary signal applied thereto with an undelayed 
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binary signal simultaneously developed by said first means 
and for developing a count representing the numerical 
difference between the number of identical and unidenti- 
cal comparisons to thereby develop the distorted auto- 
5 correlation function of said specified signal; and means 
for correcting said distorted autocorrelation function to 
develop a true autocorrelation function. 

3. The apparatus of claim 2 wherein said specified rate 
is the Nyquist rate. 

10 4 . The apparatus of claim 2 wherein said true autocor- 

relation function is equal to R xll and said distorted auto- 
correlation function is equal to 

5 sin— 1 (R xll ) 

15 2 

5. Apparatus for generating the autocorrelation func- 
tion of a specified signal comprising first means for de- 
veloping binary signals in accordance with the algebraic 
sign of said specified signal; a plurality of delay means 

20 each capable of delaying a signal applied thereto by a dif- 
ferent specified amount; sampling means for applying said 
developed binary signals to all of said delay means at a 
specified rate; a plurality of comparing means each of 
which is connected to a different unique one of said delay 
25 means and each capable of providing an output signal 
each time the delayed binary signal applied thereto is 
identical to the undelayed binary signal simultaneously 
developed by said first means; first counting means for 
counting the number of binary signals applied to all of said 
30 delay means; a plurality of second counting means each 
of which is connected to a different unique one of said 
comparing means for counting the number of output sig- 
nals provided thereby; subtraction means connected to 
said first counting means and each of said second count- 
35 ing means for determining the numerical differences be- 
between the number of identical and unidentical compari- 
sons performed by each comparison means to thereby 
develop a distorted autocorrelation function of said speci- 
fied signal; and means for correcting said distorted auto- 
40 correlation function. 

6. Apparatus for generating the autocorrelation func- 
tion of a specified signal comprising first means for de- 
veloping binary signals in accordance with the algebraic 
sign of said specified signal; a delay line 'having an input 
terminal and including a plurality of identical delay ele- 

40 ments and a plurality of taps each of which is connected 
to said line between a different unique pair of adjacent 
elements; sampling means for applying said developed 
binary signals to said delay line input terminal; a plurality 
of comparing means each of which is connected to said 
0 delay line input terminal and to a different unique one 
of said taps for developing an output signal each time the 
respective binary signals applied thereto are identical; first 
counting means for counting the number of binary signals 
r - applied to said delay line input terminal; a plurality of 
00 second counting means each of which is connected to a 
different unique one of said comparing means for count- 
ing the number of output signals developed thereby; sub- 
traction means connected to said first counting means and 
each of said second counting means for determining the 
bU numerical differences between the number of identical 
and unidentical comparisons performed by each com- 
parison means to thereby develop a distorted autocorrela- 
tion function of said specified signal; and means for cor- 
recting said distorted autocorrelation function. 

5 7. Apparatus for generating the cross correlation func- 

tion of first and second signals comprising first means for 
successively developing first binary signals, at specified 
sampling rate in accordance with the algebraic sign of said 
70 first signal; second means for successively developing sec- 
ond binary signals, at said specified sampling rate, in ac- 
cordance with the algebraic sign of said second signal; 
a plurality of third means for respectively comparing each 
of said first binary signals with a plurality of said sec- 
75 ond binary signals delay in time by different amounts for 
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determining a distorted cross correlation function of said 
first and second signals; and means for correcting said 
distorted cross correlation function. 
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